On the nonlinear dynamics of topological solitons in DNA 
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Dynamics of topological solitons describing open states in the DNA double helix are studied in the frame- 
works of the model which takes into account asymmetry of the helix. It is shown that three types of topological 
solitons can occur in the DNA double chain. Interaction between the solitons, their interactions with the chain 
inhomogeneities and stability of the solitons with respect to thermal oscillations are investigated. 
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I. INTRODUCTION 

It is widely accepted now that the DNA molecule has a 
rather moveable internal structure, and that the internal DNA 
mobility plays an important role in functioning the molecule. 
The thermal bath where the DNA molecule is usually im- 
mersed, collisions with the molecules of the solution which 
surrounds DNA, local interactions with proteins, drugs or 
some other ligands lead to activation of different types of in- 
ternal motions. Small oscillations of individual atoms near 
equilibrium positions, rotational, transverse and longitudinal 
displacements of atomic groups (phosphate groups, sugars 
and bases), motions of the double chain fragments having sev- 
eral base pairs lengths, local unwinding of the double helix, 
transitions of DNA fragments from one conformational form 
to another, for example, from A-form to B-form and so on, 
are only some of them. A more detailed list of internal mo- 
tions and of their dynamical characteristics can be found in the 
works of Fritzsche 1 1], Keepers and co-authors |2], McClure 
0, McCommon and co-authors 01, Yakushevich 00]). 

Different approaches to the modeling of the internal DNA 
mobility are known. One of them has been developed by Pro- 
hofsky and co-authors , who considered DNA as a 

lattice and took into account the motions of all atoms (except 
of hydrogen atoms) in the lattice cell. Their approach was lim- 
ited, however, by harmonic approximation, and this limitation 
did not permit them to model large amplitude internal motions 
such as, for example, local unwinding of the double helix. An- 
other approach, based on the methods of molecular dynamics 
and proposed firstly by Levitt [11] and Tidor and co-authors 
lTl2il . is known now as one of the most powerful tools of inves- 
tigation of the internal DNA mobility 1 13]. This approach is 
not limited by harmonic approximation and therefore it can be 
used to study internal motions of both large and small ampli- 
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tudes. The approach has, however, one essential deficiency: 
because of the limited possibilities of modern computers it 
can not be used to study long DNA fragments, and therefore 
it does not suitable for studies of the processes of propagation 
of local structural distortions along the molecule. 

In this paper, to investigate the internal DNA mobility, we 
use the approa ch developed in a series of works 1141 fix ua. 
EIIlIllMlinillilHIil. Peculiarity of the approach 
is that it uses rather simple models of the internal DNA dy- 
namics, which take into account only one or a few types of 
the DNA internal motions. This simplification gives an oppor- 
tunity to find analytical solutions of corresponding dynamical 
equations imitating both small and large amplitude internal 
motions. And one more merit of the approach is that it gives a 
possibility to study the internal dynamics of long DNA frag- 
ments. Three works in the series are of most interest. 

The first one has been done by Englander and co-authors 
fl4ll who studied the dynamics of DNA open states. Their 
model took into account only rotational motions of nitrous 
bases, which as it was suggested, made the main contribu- 
tion to formation of the open states. Another paper belonged 
to Peyrard and Bishop 12211 . who studied the process of DNA 
denaturation. Suggesting that the stretching of the hydrogen 
bonds in pairs made the most contribution into the process, 
they created a simplified model where only transverse mo- 
tions of bases along the direction of the hydrogen bonds where 
taken into account. The third important paper was published 
by Muto and co-authors [21]. These authors suggested that 
two types of internal motions were made the main contribu- 
tion to DNA denaturation process: transverse motions along 
the hydrogen bond direction and longitudinal motions along 
the backbone direction. Their model consisted of two polynu- 
cleotide strands linked together through the hydrogen bonds 
described by a Lennard-Jones potential, and the phosphodi- 
ester bridges in the backbone were described by an anhar- 
monic Toda potential. 

Further development of the approach was limited for sev- 
eral years by small improvements of the models and their 
combinations, and only involving numerical methods of sim- 
ulation of the internal DNA dynamics gave a new impulse and 
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FIG. 1 : Fragment of the DNA double chain consisting of 3 AT base 
pairs. Longitudinal pitch of the helix a — 3.4 A; transverse pitch 
h = 16.15 A. 



TABLE I: The values of the parameters m a , r a , I a = m^r 2 ,, for 
all possible bases a (m p — 1.67343 ■ 1CP 27 kg is the proton mass). 



a 


m a (m p ) 


r a (A) 


I a (xHT 47 m 2 kg) 


A 


135.13 


5.8 


7607.03 


T 


126.11 


4.8 


4862.28 


G 


151.14 


5.7 


8217.44 


C 


111.10 


4.7 


4106.93 



interesting possibilities which have been realized in the works 
of Van Zandt 1 26], Techera and co-authors l27ll . Salerno l25ll . 
Barbi |28, 29], and Campa 1 30] . Just these methods permitted 
not only to study a possibility of appearance of large ampli- 
tude localized distortions in the DNA structure, but also to 
investigate their stability, the influence of thermal noise, the 
interactions between the distortions, the propagation of them 
along the homogeneous and inhomogeneous DNA. 

In all these works, however, the asymmetry of the base 
pairs was neglected. That is both bases in a pair were 



H = { ^4,1^,1 + ^n,2¥£,2 + e «,i sin 2 



The first two terms of Hamiltonian correspond to the ki- 
netic energy of the n-th base pair. Here J n> i is the mo- 
ment of inertia of the n-th base of the first chain; I n jz is 
the moment of inertia of the n-th base of the second chain, 
point denotes differentiation in time t. For the base pair a(3 
(a/3 = AT, TA, CG, GC) the moment of inertia is equal 



modeled as identical structural elements with the same 
characteristics (masses, moments of inertia and so on). 
But even in the case of homogeneous (synthetic) DNA the 
asymmetry exists. Indeed, if, for example, one of the polynu- 
cleotide chains consists of only adenines, the other chain 
should consist of thymines, and this homogeneous model is 
substantially asymmetrical. Just this type of asymmetrical 
model is studied in this work. To simplify calculations, we 
consider only rotational motions of nitrous bases around 
the sugar-phosphate chains in the plane perpendicular to 
the main axis of the double chain. We find solitary wave 
solutions describing open states in the double helix. We 
classify the solitons, investigate stability of the solitons with 
respect to thermal oscillations, interactions between the 
solitons, interaction of the solitons with inhomogeneities of 
the chain. To solve all these problems, we use numerical- 
variation methods efficiency of which was proved in the 
works JH E3 13 H El Hi, devoted to the analysis of 
nonlinear dynamics of molecular chains and polymer crystals. 



II. DISCRETE MODEL OF THE DNA DOUBLE HELIX 

Let us consider B-form of the DNA molecule, the fragment 
of which is presented in Fig. ^ The lines in the figure cor- 
respond to the skeleton of the double helix, black and grey 
rectangles correspond to bases in pairs (AT and GC). Let us 
focus our attention on the rotational motions of bases around 
the sugar-phosphate chains in the plane perpendicular to the 
helix axis. Below we shall call the chain placed on the left 
by the first chain, and the right chain - by the second one. 
Positive directions of the rotations of the bases for each of the 
chains are shown in Fig. ^ 

Let us consider the plane DNA model where the chains of 
the macromolecule form two parallel straight lines placed at 
a distance h from each other, and the bases can make only 
rotation motions around their own chain, being all the time 
perpendicular to it. Let us suggest that ip rly i is the angular 
displacement of the n-th base of the first chain, and ip n> 2 is 
the angular displacement of the n-th base of the second chain. 
Then the Hamiltonian of the double chain takes the form 



Pn,l . . 2 fn+1,2 ~ <£ra,2 . T , / si 

he n , 2 sin \-V a p{ipn,l,<Pn,2)> • (1) 

I 

to I n .i = ra a r 2 a , I n £ = mpr^. The value of the base 
mass m a , the length r a and corresponding moment of iner- 
tia I a — m a r^, for all possible base pairs are presented in the 
Table|H 

The third and the fourth terms in Hamiltonian Q describe 
interaction of the neighboring bases along each of the macro- 
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molecule chains. Parameter e n> i characterizes the energy of 
interaction of the rt-th base with the (n + l)-th base of the i-th 
chain (i = 1,2). The value of the parameter is unknown. But 
if we take into account that angular displacement of one base 
is accompanied not only by overcoming the barrier due to the 
stacking interaction, but also by substantial deformation of the 
dihedral and valence angles, we can suggest that the energy of 
the displacement e n j sh ould be wittingly more than the stack- 
ing 40 60 kJ/mol 1 37], and it should weakly depend on the 
type of the base. This gives us a possibility to suggest later on 
that e„.i = e„ 2 = e > 60 kJ/mol. 

The fifth term in Hamiltonian Q corresponds to the energy 



of interaction between conjugated bases of different chains. 
Here index a(3 =AT, TA, GC, CG determines the type of the 
base pair. It is convenient to model the energy of interaction 
of conjugated pairs by the potential 



n,l ; ^n,2 ) 



l 



K a p\R. n 



R 



o |2 



(2) 



where R„ is the vector connecting the end of the base (n, 1) 
with the end of the base (n, 2), R° is the value of the vector 
for the ground state of the chain (/9„ i = 0, cp n 2 = 0. Potential 
(0 can be written in a more simple form 



J 



V a p(<pn,l,<Pn,2) = K a/3 {r a (r a + rp)(l 



cos wi) + rfj(r a + rp)(l 
I 



i<p n , 2) - r a rp[l - cos(v?„ i i - ^,2)]} 



(3) 



The rigidity of interaction K a p can be estimated from the 
energy of interaction 



e a p 



V a p(0,^) + V a p(^,0) 



1 



The pair AT (TA) is stabilized by two hydrogen bonds (they 
are shown in Fig. [flby dotted lines), and the pair CG (GC) - 
by three hydrogen bonds. Therefore we suggest later on that 
eAT = e T A = 2e C G/3 = 2e GC /3 = e. 

For the value of the energy of interaction of the bases in 
AT base pair we can take the double energy of hydrogen bond 
e = 40 kJ/mol. Then the rigidity of the bond between the 
bases is equal to 



2 2 

Kat — K T a = g-^cc = g-^CG = 



K = 



2e 



0.234 N/m. 



(4) 



'ft 



On the other hand, the value of the parameter can be estimated 
from the frequency spectrum of low amplitude oscillations of 
the chain. We shall obtain it in the next section. 



HI. DISPERSION EQUATION 

The system of equations of motion, which corresponds to 
macromolecule Hamiltonian Q, takes the form 



I<n,Wn,l — — 



dH 

dipn.l ' 

dH 



n = 0,±1,±2,. 



(5) 



Let us consider homogeneous macromolecule where only one 
type of base pairs exists a(3 (I n ,i = I a , In,2 = Ip)- 



Insert small amplitude plane wave 

(¥>n,i(*).¥>n,a(t)) = fai,¥*)4e <( ' OB - wl) , 

into the system of equations (|5}. Here (ipx, if 2) is a normal- 
ized to 1 two dimension vector, A <C 7r is the amplitude, 
q 6 [0,%/a] is the wave number. It is easy to show that in 
the linear approximation the frequency u> should satisfy the 
dispersion equation 



Boj 2 + C = 



(6) 



where 



B = [K a p{I a rj + Iprl) + 4«(/ Q + Ip) sin 2 ^)}/I a Ip ; 
C = [AnK a0 {rl + r%) sin 2 ^ + 16« 2 sin 4 ^]/I a Ip , 

k = e/2 is the rigidity of the interaction of neighboring bases 
along the chain. 

Dispersion curve (|5Jl has two branches: 



[{B- VS 2 -4C)/2] 1 / 2 ; 
[(B + V / S 2 -4C)/2] 1 / 2 . 



The upper curve co = Lo {q) corresponds to optical phonons, 
the lower curve u> — u> a (q) corresponds to acoustic phonons 
in the chain. 

The frequency u> a (q) tends to zero as q — > 0. Let us deter- 
mine the velocity of acoustic phonons as 



i-'O 



Urn 

q^0 q 



I a rl + Iprl 



The dependence of the sound velocity vq in the homogeneous 
molecule a(3 =AT (GC) on the energy of rotation e is pre- 
sented in the TablelTTl 

According to different estimations 1 38, 39, 40] the veloc- 
ity of sound in DNA is on the interval from 1890 m/s till 







71/2 7C 

qa 



n/2 7C 
qa 



FIG. 2: Acoustic to — Lu a (q) (curve 1, 3, 5) and optical u = tj (q) (curve 2, 4, 6) branches of the dispersion curve for homogeneous chain 
(a/3 =AT, e =60, 600, 6000 kj/mol) for chain with K = 0.234 N/m (a) and K = 4.744 N/m (b). 



3500 m/s. From the Table |ll] it is clear, that among three 
typical values e = 60, 600, 6000 kJ/mol the value e = 
6000 kJ/mol is the best one. Just this value will be used in 
the numerical investigations of the dynamics of topological 
soli tons. 

The lowest value of the optical frequency is 

w o (0) = yjK a0 (I a r} + I rl)/I a I p . (7) 



TABLE II: Dependence of sound velocity Vo (m/s) on the value of 
the parameter e for homogeneous a/3 =AT (GC) chain 



e (kJ/mol) 


60 


600 


6000 


AT 


219.47 


694.02 


2194.7 


GC 


223.38 


706.39 


2233.4 



According to J4lll w o (0) = 35 cm 1 , therefore from we 
have 

Kat = K = 4.744 N/m, K C g = \ K = 7 - 117 N/m - ( 8 ) 

This estimation of the value of the rigidity differs from that 
obtained in (|4}. The use of the value given in (0} gives sub- 
stantially lower value of the frequency u> o (0) = 7.77 cm -1 . 
Thus we have the following estimation of the value of the pa- 
rameter A": 0.234 < K < 4.744 N/m. The view of the disper- 
sion curves for homogeneous chain (afi = AT) with different 
values of the parameters is presented in Fig. |2] 

For numerical investigation of the soliton dynamics we 
shall take intermediate value K = 0.8714 N/m which cor- 
responds to the frequency w o (0) = 15 cm -1 , and energy of 
interaction eat = 149 kJ/mol. 
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IV. NUMERICAL METHOD OF FINDING SOLITARY 
WAVE SOLUTIONS 

Complexity of the system of equations of motions (0 does 
not permit us to carry out analytical investigation. Therefore, 
we shall study it numerically and use variation technique, pro- 
posed in 13211 - to find soli ton like solutions. 

Let us consider homogeneous DNA molecule (for all n 
In,i = la, I n ,2 = Ip, where af3 = AT (TA, CG, GQ). We 
shall find the solution of system Q in the form of a wave with 
smooth constant profile. For the purpose, let us suggest that 
l Pn,i(t) = ^i(£)> fn,2{t) = <P2{0> where the wave variable 
£ = na — vt, and v is the velocity of the wave. 

Let us assume, that the functions ip\ and ip 2 smoothly de- 
pend on £. Then the time second derivatives can be substituted 
for discreet derivatives 

(P(f„,i 2 dlfii 2 2 

~ V ■^To = V VPn+\,i ~ 2<fn,i + (Pn-l,i)/a , (9) 



dt 2 



i = 1,2. Using these relations, we can write the equations of 
motions (0 in the form 



L ViiX =Q, L PB , a =0, n = 0,±l,±2,. 
Here the functional 



(10) 



v 2 



-e(sm '— +sm — -) 

-V a p(<p n.l i ¥>n,2)} 



is a discreet version of the Lagrangian 



-e(sin 



n ~ 

■ 2 fn+1,1 — <Pn,l ■ 2 fn+1,2 ~ </?ra,2 



-V a /3{ip n ,l,<Pn,2)] , 



which corresponds to the system of equations of motion 0. 

For further analysis it is convenient to write the functional 
L in the dimensionless form 

L = 2L/K(r 2 A + r 2 T ) = 

/\c a (<p n +l,l ~ Pn,l) 2 + Cp(tp n +1,2 ~ Vn^) 2 
n 

2 fn+1,1 - <Pn,l , ■ 2 fn+1,2 ~ <Pn,2 ■, 

-g{sm +sin ) (11) 

-U a p(ip nt u <Pn,2)], 

where the dimensionless coefficients 

v 2 I a v 2 I a 



Ka 2 {r 2 A +r 2 T Y ° Ka 2 {r\ +r^) : 



parameter of cooperativity 

g = 2e/K(r 2 A + r 2 ), (12) 
dimensionless potential U a p{f n ^,tp n .2) = 
2V a p(<p nt i,ip nt 2)/K(r A + r|). 

Soliton solution of the system fl!0l > can be found numeri- 
cally as a solution of the problem on conditional minimum 



mm 

<P2,i,---,<PN-l,i, 1=1,2 



<Pl,l = f-oo,l, <Pl,2 = f-oo,2i 
<PN,1 = foo,!, <pN,2 = ¥>oo,2 



(13) 

(14) 
(15) 



Boundary conditions 1141 . dl5> for the problem Jl 3I > deter- 
mine the type of the soliton solution. We should take a rather 
large number N, in order that the form of the solution of the 
problem might not depend on its value. For the purpose it is 
enough to take N ten times larger than the width of the soliton. 

The soliton solution of the problem Jl 31 can be character- 
ized by the topological charge q=(<2i, 92), where = (foo,i — 
<P-oo,i)/27r, i — 1, 2, is an integer (gi = 0, ±1, ±2, ...). To 
find soliton solution with topological charge q, it is necessary 
to solve the problem on minimum M3\ with boundary condi- 
tions 

^-00,1 = ^-00,2 = 0, <A»,i = 27rgi, ^00,2 = 27rq 2 ■ 

This problem was solved by the method of conjugated gradi- 
ent. The value N = 2000 was taken, and the initial point 

<Pn,i = [1 + tarda fi(n - N/2)]nq l , i = 1,2 , 

was used. Here fj, is a changeable parameter. 

Soliton solution (solution in the form of a solitary wave) 
{<^° x>9f° 2}n=i corresponds to topological soliton with the 
energy E = K(r\ + r\~)Ej2, where dimensionless energy 

N-l 

E = ^ [Ca&n+l,! - <Pns) 2 + Cp(<Pn+l,2 ~ ^n^) 2 + 



n=l 



. / . 2 Vn+1,1 - <Pn,l . ■ 2 fn+1,2 ~ ¥n,2 , 
+<7(sm ■ h sin ■ — 



+ U a p(ip ny iip n ^)} 



and with the diameter 



where the point 



D = l + 2^{n-n) 2 Pn , 



N 

71=1 



determines the position of the soliton center, and the formula 
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p„ = E n /E = <j i [c a ((f n+ltl - fn-l,l) 2 + Cp((p n+ x,2 ~ fn-1,2) 2 ] + 7^9 



■ 1 fn+1,2 — fn,2 , . 2 Vn.l ~ Vn— 1.1 , • 2 ^",2 — Vn-1,2 

sm ; ■ — h sin — ■ ■ — h sin — ■ — 



. 2 fn+1,1 — fn.l 

sin : — 



+ U a /3{lfn,lVn,2) > /E 



r 



gives the distribution of the energy along the chain. 



V. DYNAMICAL PROPERTIES OF SOLITONS 

At the beginning let us consider stationary soliton solutions 
of the problem Jl 31 . In the dimensionless functional L, coeffi- 
cients c a = eg — when v = 0. So, only one dimensionless 
parameter g d!2t which characterizes cooperativity of rota- 
tional motions remains in functional il It . The existence of 
soliton solution and its form depend on the value of the pa- 
rameter. 



A. Stationary solution 

The results of numerical investigations of the problem Jl 3i 
show that in the homogeneous chains stationary topological 
soliton solutions exist when the parameter of cooperativity g 
is larger than the threshold value: g > go > 0. The ab- 
sence of the soliton topological stability when g < go can 
be explained in the following way. Any topological defect 
can be eliminated by turning the DNA bases. The turning of 
one base about 360 degrees transfers the system to the initial 
state, that is <fri y „ = 0j jn ± 2tt. And this is why the narrow 
solitons with the size equal to one link of the chain are equiv- 
alent to the ground state and this is why they are unstable. So, 
only relatively wide solitons with the general turning consist- 
ing of several small changes of rotational angles that is when 
\4>i,n+l ~ 4>i.n I ^ 27T, are stable. Dependence of the threshold 
value go on the soliton topological charge q for homogeneous 
AT and GC chains are given in the Table lllll 

From Fig. [3] it becomes clear that the soliton energy E 
and its width D monotonously increase when the parameter 
of cooperativity g increases. For soliton stability it is neces- 
sary that its width D > 4.33. The view of stationary solitons 
with the parameters of cooperativity g = 10 and g = 150 
is presented in Fig. |4] In the case of the soliton with topo- 
logical charge q=(l,0), the first component has the form of a 
smooth step (when n monotonously changes, the base of the 



TABLE III: Dependence of the threshold value of the parameter of 
cooperativity go on the value of the soliton topological charge q= 
(?i, 12) for homogeneous a/3 =AT (GC) chain 



first of two DNA chains makes a complete turn) accompanied 
by smooth small amplitude deformation in the second compo- 
nent (Fig. @^). In the case of the soliton with q=(0,l), only 
the second component has the form of a step (Fig. |4j>). In the 
case of soliton with q=(l,l), each of the components has the 
form of a step (Fig. gj;), steps being displaced relatively one 
another. Later on we shall show that this soliton is the bound 
state of two topological solitons with the charges qi=(l,0) and 
q2=(0,l). There exist two equivalent states of the soliton: the 
left state q=(l,l);, when the soliton with the charge qi is on 
left of the soliton with the charge q2 (Fig. and the right 
state q=(l,l) r , when the soliton with the charge q! is on right 
of the soliton with the charge q2. 

When e = 6000 kJ/mol and K = 0.234 N/m the pa- 
rameter of cooperativity g = 150.24 ^> go, and when 
K = 0.8714 N/m the parameter g = 40.34 > g (see Ta- 
ble lllH for all types of topological solitons. So, for these val- 
ues of the rigidity parameter K stable solitons with different 
topological charges exist. For maximum value of the rigid- 
ity parameter K = 4.744 N/m the parameter of cooperativity 
g = 7.41 < go, and this means that stable topological solitons 
are absent. Thus the problem of the existence of topological 
solitons (open states) in the DNA reduces to the problem of 
receiving exact estimation of the parameters e K. But this is, 
however, a rather difficult problem. We think that the values 
e = 6000 kJ/mol, K = 0.8714 N/m, when all three types of 
solitons exist, are the most grounded. And we shall use these 
values for further calculations. 

Dependence of the energy E and the diameter D of sta- 
tionary topological soliton on its topological charge q in the 
chain with e = 6000 kJ/mol and the transverse rigidity de- 
termined by formulas @ and (|8j, are given in the Table IIVI 
From the data of the table it follows that the energy of the 
interaction of the solitons with the charges (1,0) and (0,1) is 
equal to AE = E(1,0) + E(0, 1) - E(l, 1) = 90.71 kJ/mol 



TABLE IV: Dependence of the energy E and the diameter D of 
stationary topological soliton on its topological charge q with two 
values of the transverse rigidity K 



(gl,g2) 

AT 
GC 



(LP) 
8.3 
12.0 



(0,1) 
5.7 
8.2 



(1,1) 
8.3 
12.0 







af3 = 


AT 


a/3 = 


GC 


K (N/m) 


q 


E (kJ/mol) 


D 


E (kJ/mol) 


D 




(1,0) 


2776.52 


16.59 


3405.48 


13.69 


0.234 


(0,1) 


2237.36 


19.58 


2733.60 


16.18 




(1,1) 


4971.72 


42.85 


6087.54 


35.13 




(1,0) 


5329.76 


9.03 


6394.96 


7.64 


0.8714 


(0,1) 


4302.09 


10.59 


5146.38 


8.98 




(1,1) 


9551.14 


22.60 


1 1444.96 


18.93 



150 



100 



icq 



50 



(a) 


^3 







150 



100 





FIG. 3: Dependence of the dimensionless soliton energy E and the width D on the velue of the parameter of cooperativity g when q=(l,0), 
(0,1), (1,1) (curves 1, 2, 3) in homogeneous AT (a),(c) and GC chains (b),(d). 



for the chain with the rigidity of the transverse interaction 
K = 0.8714 N/m, and the energy is equal to AE = 42.16 
kJ/mol for the chain with K = 0.234 N/m. 



B. Nonstationary solutions 

Numerical investigation of the problem Jl 3I > shows, that in 
the homogeneous chain, topological soliton has the interval of 
the velocities < s < s\ < 1, where s = v/vq is the dimen- 
sionless velocity and vo is the velocity of sound. Dependence 
of maximum velocity of the soliton si on its charge q, on type 
of the base in the chain aj3 and on the rigidity of the transverse 
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FIG. 4: The view of stationary soliton with the topological charge 
q=(l,0) (a); q=(0,l) (b); q=(l,l); (c). Continuous lines correspond 
to displacements by the first component y> n ,\ \ dotted lines - to dis- 
placements by the second component (f n ,2\ thin lines correspond to 
chain with g = 10; flat lines - to chain with g — 150. 
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FIG. 5: Dependence of the energy E and the diameter D of the soli- 
ton which moves along homogeneous AT chain and has the charge 
q= (1, 0) (curves 1 and 4), q= (0, 1) (curves 2 and 5) and q= (1, 1) 
(curves 3 and 6), on the dimensionless velocity s (e = 6000 kj/mol, 
K = 0.8714 N/m). 



TABLE V: Dependence of the maximum value of the soliton veloc- 
ity si on its topological charge q, the soliton moving in the homoge- 
neous a/3 chain with the transverse rigidity K. 



K (N/m) 


q 


AT 


GC 




(1,0) 


0.77 


0.70 


0.234 


(0,1) 


0.88 


0.84 




(i,D 


0.77 


0.70 




(1,0) 


0.64 


0.55 


0.8714 


(0,1) 


0.86 


0.84 




(1,1) 


0.65 


0.56 



TABLE VI: The dependence of the soliton mass of rest M (the val- 
ues are given in proton mass units m p ) on its topological charge q, 
the soliton being in the homogeneous a/3 chain with the transverse 
rigidity K. 



K (N/m) 


q 


AT 


GC 




(1,0) 


7640 


9663 


0.234 


(0,1) 


4052 


4064 




(i,D 


11581 


13590 




(1,0) 


14978 


18730 


0.8714 


(0,1) 


7899 


7804 




(1,1) 


22660 


26265 



interaction K is presented in the Table IV! 

Dependence of the soliton energy E and the diameter D on 
the dimensionless velocity s is presented in Fig. [5] With the 
increasing of the soliton velocity its energy monotonically in- 
creases, and the diameter monotonically decreases. Using the 
dependence E(s) we can find the mass of rest of the topolog- 
ical soliton 

M = Um W£)_l«. 

Dependence of the mass of rest M of the soliton on its charge 



q, on the type of the bases of the chain a/3 and on the rigidity 
of the transverse interaction K is presented in the Table lVIl 

Numerical investigation shows that all topological solitons 
at all permitted velocities are stable. They move along the 
chain with constant velocity, their form and energy being 
conserved. Thus, the specificity of the chain of the DNA 
molecule leads to a principal effect consisting in the possi- 
bility of preferable localization of soliton excitations on one 
chain. Moreover, it appears that soliton excitation with the 
charge (1,1) is a bound state of two excitations localized in 
separate chains. 
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FIG. 6: Potential of interaction of solitons C/ qi , q2 (R) with qi = 
(1, 0), q 2 = (0, 1) (curve 1); qi = (1, 0), q 2 = (0, -1) (curve 2); 
qi = q2 = (1, 0) and qi = q 2 = (0, 1) (curve 3). 

VI. INTERACTION OF TOPOLOGICAL SOLITONS 

DNA is a rather long molecule, and several open states can 
be activated in it simultaneously. Therefore it is interesting to 
consider the problem of interaction of solitons imitating the 
open states. 

Numerical approach d 1 3I > permits us to investigate the prob- 
lem and to obtain the dependence of the energy of a pair of 
solitons with the charges qi, q2 on the distance between then- 
centers rii, ri2- For the purpose, it is necessary to take the 
boundary conditions and the initial point, which correspond 
to a pair of topological solitons with the centers moving away 
at a distance R. When minimizing the energy of the system 
E = — L, it is necessary also to fix the turns of the bases, 
which correspond to the centers of solitons. Then the energy 
of the obtained state E(R) corresponds to the energy of a 
pair of topological solitons, the solitons being at a distance 
R = n% — n\ from each other. By changing positions of the 
soliton centers, we can obtain the potential of interaction 

U nuCl2 (R) = E(R) - E( qi ) - E(q 2 ) , 

where E(qx) and E(q%) is the energy of isolated solitons. 

The potential of interaction of solitons of different types 
and the potential of topological solitons with the charges of 
the same sign are presented in Fig. [6] The potential of two 
solitons of different types with qi = (1, 0) and q-2 — (0, 1) 
has the form of symmetrical double well potential (Fig. |6j 
curve 1). Maximum of the potential is reached when R = 0, 
that is when the centers of the solitons are placed at the neigh- 
boring chains and when the solitons are opposite to each other. 
From energetic point of view this configuration of the solitons 
of different chains is the most disadvantageous. Minimum of 
the energy is reached when R — ±20a. Thus two solitons of 
this type can form two energetically equivalent coupled states. 
One of the states (left minimum of the potential of interaction) 
corresponds to the left isomer of the topological soliton with 
the charge q= (1, l)z, and the other state (the right minimum 
of the potential of interaction) corresponds to the right isomer 
of the soliton q = (1. l) r . 

If solitons have different signs of charges qi = (1,0) and 



q2 = (0, —1), the potential of interaction has a bell-like form 
with one maximum at R = (Fig. |5J curve 2). From the 
potential it follows that, the solitons which belong to different 
chains, should repulse from each other. Solitons with the same 
sign of charge (qi = q2 = (1,0), (0,1)) also repulse from 
each other. 

When distance between the solitons decreases, the energy 
monotonically increases and goes to infinity when R — > 
(Fig.|6j curve 3). 

The potential of interaction C/ qi iqa (R) permits to predict 
the result of repulsion of solitons with the charges qi and 
q2. Let us model the repulsion of the solitons. For the pur- 
pose, let us consider a double chain consisting of N = 4000 
base pairs. At the ends of each of the polynucleotide chains 
let us introduce viscous friction which provides with absorp- 
tion of phonons. The system of equations of motion 10, 
n = 1,2, ...,7V, was integrated numerically with the initial 
condition which corresponds to two topological solitons with 
the centers placed in the points n\ = N/A and ri2 = 37V/4 
and with the velocities si = — s 2 = s > 0. 

The results show that collision of solitons having equal 
signs qi = q2 = (±1,0), (0, ±1) leads to their reflection 
at one another. When the velocities s are small, the reflection 
is practically elastic, and when the velocities s are large, colli- 
sion is accompanied by slight emission of phonons. Collision 
of solitons with qi = (±1,0), q2 = (0, =Fl) and s = 0.5 
leads to their reflection, accompanied by slight emission of 
phonons. This behavior is in a good agreement with the form 
of corresponding potential of interaction (Fig. [6] curve 2). To 
pass through one another, solitons need to overcome energy 
barrier C/ qii q 2 (0) = 2025 kJ/mol. Thus, their kinetic energy 
should be equal to Ek(s) = Eq 1 (s) + Eq 2 (s) — S qi (0) — 
Eq 2 (0) > f7 qi , q2 (0). This condition is fulfilled only in the 
vicinity of the most possible values of the velocity (see. Fig. 
[5}- So, when s = 0.5, the kinetic energy Ek — 1732 kJ/mol 
is lower than the height of the energy barrier (reflection takes 
place), and when s = 0.6, the energy Ek — 3000.9 kJ/mol is 
higher than the barrier (solitons pass through one another). 

Solitons with the charges qi = (±1,0), q2 = (0, ±1) at- 
tract one another at a distance R > 20a, and when the dis- 
tance is shorter they repulse one another. Here the energy 
barrier [7 qi iq2 (0) = 4989 kJ/mol does not permit solitons to 
pass through one another. Solitons always reflect. Formation 
of the bound state does not occur even when the value of the 
velocity is small. It is explained by small value of the bond 
energy AE — 91 kJ/mol. 

Potential of interaction of different topological solitons 
with the charges of opposite signs (qi = — q2) is presented 
in Fig. Potential of interaction of one component solitons 
(qi = (1,0), (0, 1)) monotonically decreases with the de- 
creasing of the distance between the solitons. When R — * 0, 
the potential J7 qi , q2 (i?) — » —[E(qi) + i?(q 2 )]. At a distance 
R = solitons completely recombine. 

Potential of interaction of two component solitons with the 
charges of different signs has a similar form. If the solitons 
have different polarity, that is, if the first soliton is a left isomer 
(qi = (1, 1)/) and the second soliton is a right isomer (q2 = 
(—1,-1), ). Solitons of that type attract one another. Their 
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FIG. 7: Potential of interaction of solitons Uq 1 ,q 2 (R) with qi = 
(1, q 2 = (-1, -1)/ (curve 1); qi = (1, 1),, q 2 = (-1, -l) r 
(curve 2); q x = (1, 0), q 2 = (-1, 0) (curve 3); qi = (0, 1), q 2 = 
(0,-1) (curve 4). 



collisions always leads to recombination of the solitons. 

If solitons have the same polarity, they repulse when R > 
20a, and attract when the distances are shorter. Recombi- 
nation of the solitons requires overcoming the energy bar- 
rier 1730 kJ/mol. Solitons overcome the barrier only when 
s > 0.38, and the value of their kinetic energy is more than the 
height of the barrier. When the value of the velocity is smaller, 
solitons reflect, and when the value is larger they recombine 
(Fig.|§J. During recombination the energy of solitons is spent 
for intensive emission of phonons. Breather-like excitations 
can be also formed. 

Collision of one component soliton with two component 
soliton can lead (depending on the relationship of the signs 
of charges and on the polarity of the two component soliton) 
to their partial recombination or to inelastic reflection accom- 
panied by disintegration of the two component soliton. So, 
at the velocity s = 0.5 collision of the soliton having the 
charge qi = (1,0), with the two component soliton having 
the charge q2 = (1,1)/, leads to inelastic reflection of the first 
soliton, and at the same time the second soliton disintegrates 
into two one component solitons with the charges (1,0) and 
(0,1) (Fig. 0. 



VII. EFFECT OF THE CHAIN INHOMOGENEITIES ON 
THE DYNAMICS OF TOPOLOGICAL SOLITONS 

Till now our investigation was limited by consideration of 
homogeneous model of DNA. The real DNA is, however, a 
substantially inhomogeneous system, therefore it is of special 
interest to consider the effect of the chain inhomogeneity on 
the dynamics of topological solitons. In the inhomogeneous 
chain, the energy of stationary soliton E will depend on the 
position of the center of the soliton n. To move, soliton re- 
quires to overcome the energetic potential barrier E(n). To 
find the energy of the soliton with the center at the point n = h 
we need to solve numerically the problem of minimizing 

E -> min : (16) 

C^2,i,...,<^JV — 1,2 






FIG. 8: Recombination of solitons with the same signs of charges 
and the same polarities (qi = (1, 1);, q 2 = ( — 1, —1);, s = 0.5). 

<fil,l = f-oo,l> <Pl,2 = ¥>-oo,2, (17) 
<PN,1 = <A»,1, <^AT,2 = </2oo,2, (18) 

where the energy 
E = 2J e ( sm 2 2 ' 

n 

n,l ) ^Pn,2 )\ ) 

a(3 n - the sequence of the base pairs along the chain. The 
boundary conditions d!7i . (II 81 are the same as those in the 
problem (I13> . To fix the soliton center position we need to 
solve the problem of minimizing d!6i with respect to variables 
fn,i where 2 < n < N — 1, n ^ n, i = 1, 2. For soliton with 
charge qi ^ it is necessary to fix the value (p n j = irqi, and 
for soliton with q2 ^ - the value ipn t 2 = nq2- 

The problem on conditional minimum d!6l > has been solved 
by the numerical method of conjugate gradient. We took N = 
2000. Possible view of energetic profile of the soliton moving 
in the inhomogeneous chain is presented in Fig. ^| 

At the beginning, let us estimate the effect of point inhomo- 
geneities. From Fig. HOh it is obvious that one point defect in 
the homogeneous AT chain leads to appearance of localized 
potential barrier with height equal to Ed — 150 kJ/mol. To 
overcome the barrier, the soliton kinetic energy should satisfy 
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FIG. 9: Disintegration of two component soliton (q2 = (1,1);) 
when it comes into collision with one component soliton (qi = 
(1,0)). The velocity of the movement is s = 0.5. 



the condition Ek(s) — E(s) — E(0) > Ed- Soliton can to 
overcome the barrier only when its velocity s > Sd, where 
the threshold value of the velocity Sd is taken from equation 
E(sd) — E(0) = Ed- From the data of Fig. [5] it is easy to 
find, that for soliton with q=(l,0) the velocity Sd — 0.21, for 
soliton with q=(0,l) the velocity Sd = 0.28, and when q=( 1,1) 
the velocity Sd = 0.17. 

Let us model numerically the interaction of soliton with lo- 
cal defect of the homogeneous AT chain. For the purpose, let 
us consider homogeneous AT chain consisting of A = 4000 
bases with one base GC in the middle of the chain in the point 
n = N/2. Suggest that at the initial time a topological soliton 
is in the point n — N/ 4 and consider its movement through 
the chain inhomogeneity. The results of numerical modeling 
of the soliton dynamics show, that independently on the value 
of topological charge q soliton with velocity s = 0.05 reflect 
from this point defect, but for s = 0.5, moves through the 
point defect with negligibly small energy loss. 

Point defect in the homogeneous GC chain leads to forma- 
tion of localized potential well with depth 150 kJ/mol (see 
Fig. 110b). Almost at all values of the velocity soliton easily 
propagates along this chain without formation of bound state. 
Thus we can conclude that soliton moving in the DNA chain 
with sufficiently large velocity (s > Sd) is stable with respect 
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FIG. 10: The view of energetic relief for soliton in the inhomoge- 
neous chain: in homogeneous AT chain with one GC base pair (a); 
in homogeneous GC chain with one AT base pair (b); in the chain 
the first part of which consists of only AT base pairs, and the second 
- of only GC base pairs (c); in the chain with random sequence of 
base pairs (d). Dotted line shows the relief of the soliton with topo- 
logical charge q= (1,0), firm line - the relief for soliton with charge 

q=(M). 



to point defects. 

In the chain one part of which consists of only AT base 
pairs, and the other - of only GC base pairs energetic bar- 
rier takes the form of smooth step (Fig. UOt). The height of 
the step is equal to the difference between the values of soli- 
ton energy in homogeneous GC and AT chains. From data 
of table IIVI it follows, that the height of the step is equal 
to AE = 1065 kJ/mol for soliton with topological charge 
q=(l,0) and AE = 844 kJ/mol for soliton with q=(0,l), and 
when q=(l,l) the energy A_E=1894 kJ/mol. Soliton moving 
along homogeneous AT region of the chain, can enter into 
GC region only if its kinetic energy E k (s) > AE. As seen 
from Fig. |5] this condition is satisfied only if soliton veloc- 
ity s > Sk, where the threshold value of the velocity is de- 
termined by the equation Ek(sk) — AE. For soliton with 
q=(l,0) the velocity s k = 0.48, for q=(0,l) s k = 0.59, and 
for q=(l,l) the velocity s k — 0.52. 

Let us model numerically the soliton moving from homo- 
geneous AT region of the chain to homogeneous GC region. 
The results of the modeling show that the soliton with veloc- 
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FIG. 11: Movement of soliton with q= (1,0), s = 0.5 through the 
boundary between homogeneous AT and GC regions. 
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FIG. 12: Entering two component soliton (charge q= (1, 1);, veloc- 
ity s = 0.5) the random inhomogeneous region of the chain, and 
further disintegration of the soliton. 



ity s — 0.05 < Sfc and with any topological charge reflects 
elastically from the boundary between the regions. At given 
velocity the soliton kinetic energy is not large enough to over- 
come energetic barrier (Ek(s) <C AE). Soliton with s = 0.5 
q = (1,0) moves through the boundary between homoge- 
neous regions and its motion is accompanied by emission of 
phonons. Inside the region consisting of GC base pairs, soli- 
ton continues to move, but with a smaller magnitude of the 
velocity (Fig. II It . For given value of q the threshold value of 
the velocity = 0.48 < 0.5. So, the kinetic energy of soli- 
ton is large enough to overcome energetic barrier. Because 
the main part of the kinetic energy is spent to overcome the 
barrier, the velocity of the soliton substantially decreases after 
overcoming the barrier. When q=(0,l) the threshold value of 
the velocity Sk — 0.59 and soliton reflects at s = 0.5 from the 
boundary of the homogeneous regions. The reflection is ac- 
companied by phonon emission. For soliton with the charge 
q=(l,l) the velocity s = 0.5 < = 0.52 is not enough to 
overcome energetic barrier. Collision of soliton having topo- 
logical charge (q = (1,1);) with the boundary between the 
homogeneous regions leads to the disintegration of the soli- 
ton. It disintegrates into two one component solitons with the 
charges qi = (1,0) and q2 = (0, 1). Soliton with the charge 
q2 continues to move into GC region of the chain, and soliton 



with the charge qi reflects from the boundary. 

Let us consider the propagation of soliton in the inhomo- 
geneous chain with random sequence of bases. In this case 
random energetic relief E(n) is formed. The amplitude of the 
relief for soliton with q=(l,0) reaches 1000 kJ/mol, and for 
soliton with q=(l,l) - 1500 kJ/mol (Fig. [Tohl. It is obvious 
that uniform propagation of soliton in the chain of that type 
is impossible, because soliton loses part of energy for phonon 
emission when crossing each homogeneity. 

Let us consider the movement of soliton through inhomo- 
geneous region of the chain. For the purpose, let us suggest 
that the second part of the chain is formed by a random equal- 
possible sequence of base pairs AT, TA, CG, GC. The results 
of numerical modeling of the soliton dynamics show that soli- 
ton with small value of the velocity s = 0.05 and with any 
topological charge reflects from the boundary of the inhomo- 
geneous region. This points out, that penetration of the soli- 
ton into the inhomogeneous region requires the overcoming 
of some energy barrier. Soliton with larger velocity s = 0.5 
and charge q= (1, 0) overcomes this barrier, enters the disor- 
dered region of the chain and stops there. The movement in 
the disordered region is accompanied by intensive emission of 
phonons, which leads to the stop of the soliton. Soliton with 
q=(0, 1) can not overcome the barrier even at this value of the 
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velocity. The soliton reflects from the boundary of the inho- 
mogeneous region. The reflection is accompanied by emis- 
sion of phonons. Two component soliton with q= (1,1)/ 
enters inhomogeneous region, and at the same time it dis- 
integrates into two one component solitons with the charges 
qi = (1, 0) and q2 = (0, 1). The solitons moves some time in 
the inhomogeneous chain, then they stop (Fig. I12l> . The path 
of the solitons can reach several hundred base pairs. 

Analogous results have been obtained even in the case when 
inhomogeneous region was formed by the base pairs AT and 
TA. Thus, the sequence of nitrous bases of DNA molecule 
should substantially influence the characteristics of the motion 
of topological soliton. Note, that it has been pointed out firstly 
in the work 1 231 • 




VIII. INTERACTION OF TOPOLOGICAL SOLITONS 
WITH THERMAL OSCILLATIONS OF THE CHAIN 

Dynamics of a thermalized chain consisting of N sites, is 
described by the system of the Langevin equations 



dH 
dH 

' 2 dfn,2 

n = 1, 2, ...N 



Iri,\(pn.\ — — 



£,n,l — r/ nj i(/9 n! i ; 



in,2 - TI nt 2lpn,2 , (19) 



where the Hamiltonian of the system H is given by Eq. Q, 
£ nj i are random normally distributed forces describing the in- 
teraction of the n-th base of the i-th chain (i = 1,2) with 
thermal bath, V = l/t r is the coefficient of friction, t r being 
the relaxation time of the rotation velocity of one base. The 
random forces ^ have normal distribution and the correla- 
tion functions are 



(£n,i(*l)£m,j(*2)) = 2Tk B TS nm S ij S(t 1 - t 2 )^/ln,Jn,j, 

n,m = 1,2, ...,N, i,j = 1,2 , 

where ks is Boltzmann's constant and T is the temperature of 
thermal bath. 

The system d!9l > was integrated numerically by the standard 
fourth-order Runge-Kutta method with constant step of inte- 
gration At. The delta function was represented as S(t) = 
when|t| > At/2, and 8(t) = 1/Atwhen|t| < At/2, i.e., the 
step of numerical integration corresponded to the correlation 
time of the random force. In order to use the Langevin equa- 
tion, it was necessary to suggest that At <C t r . Therefore we 
chose At = 0.001 ps and the relaxation time t r > 1 ps. 

Let us check stability of topological soliton with respect 
to thermal oscillations of the chain. For the purpose, let 
us consider homogeneous periodical AT chain consisting of 
N = 4000 base pairs at the temperature T = 300K. Let us 
integrate system dl 91 with the initial condition corresponding 
to topological soliton (s = 0.5) with center placed in the point 
n = N/4. Numerical integration shows stability of solitons 
at all values of the charge and at both values of the transverse 
rigidity K = 0.234 N/m and K = 0.8714 N/m. The viscos- 
ity of the environment leads to quick stop of the soliton, and 
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FIG. 13: Stability of topological soliton (q= (1, 1);, s = 0.5) in 
the thermalized homogeneous AT chain (T — 300K, t r — 1 ps). 
Dependence of the distribution of angular displacements <p n ,i, tfn,2 
and energy E n along the chain on time t is shown (K = 0.234 N/m). 



after that all time it remains immovable. Soliton remains sta- 
ble with respect to thermal oscillations during the all time of 
numerical integration t = 5 x 10 3 ps ("Fig. 1131. 

Let us note that in contrast to the models of phi-4 and of 
sine-Gordon the stability of solitons in the DNA model has not 
topological nature. Solitons can be destroyed. To show this, it 
is enough to suggest that the soliton width is equal to one base 
pair (soliton of that type is equivalent to the ground state of the 
chain). Here the stability is associated with energetic factors. 
From Fig. ^]it is well seen that in the region of localization of 
the soliton, the density of the energy is equal to E n > fc^T. 

Soliton path length in the thermalized homogeneous chain 
(T = 300K) depends on the value of relaxation time t r (on 
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FIG. 14: The braking of topological soliton (charge q= (1, 1);, ini- 
tial velocity s = 0.5) in the thermalized cyclic homogeneous AT 
chain (T = 300K, K = 0.234 N/m and t r = lOOOps). 



the viscosity of the surrounding of the molecule). At strong 
viscosity t r = 1 ps soliton has time to pass only 7 chain links 
till full stop. Then it remains immovable all the time (Fig. II 3I >. 
When the viscosity is lower t r = 10 ps soliton has time to pass 
41 links, and when t r = 100 ps - 480 links. The braking of 
soliton at low viscosity (t r = 1000 ps) is shown in Fig. ^] 
Soliton passes more than 3000 chain links, and then it begins 
to move as a massive Brownian particle. 

The braking of soliton in the homogeneous chain is con- 
ditioned only by viscosity. When the viscosity is absent 
(i r = oo) soliton is moving along thermalized chain with 
constant velocity (Fig. U5t.Thetmal phonons by themselves 
do not influence the soliton dynamics. 

Let us note that topological soliton can move along the 
DNA chain in the presence of viscosity too. To organize the 
propagation it is necessary to select in a special way the se- 
quence of bases. If concentration of AT base pairs monotoni- 
cally increases, inclined potential E(n) is formed. The energy 
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FIG. 15: Movement of topological soliton (charge q= (1,1);, initial 
velocity s = 0.5) in the thermalized cyclic homogeneous AT chain 
(T = 300K, K = 0.234 N/m and t r = oo). 



overfall can reach 1 1 16 kJ/mol at K = 0.234 N/m and 1894 
kJ/mol at K = 0.8714 N/m. Soliton will propagate along the 
relief inclination as a Brownian particle moving in the viscous 
media under the action of external constant force. 

Thermal phonons substantially influence the interaction of 
soli tons. In the work 1 42], it was shown that topological soli- 
tons of the model 0-4, can interact with one another through 
thermal phonons. This interaction comes to repulsion of the 
soli tons. As a result, in the thermalized chain the interaction 
of the solitons of different charges substantially changes. At 
a long distance they will repulse. To model this phenomenon, 
let us consider collision of solitons with different charges and 
polarities in the thermalized cyclic AT chain (qi = (1,1);, 
q2 = (— 1, — l) r , si = —S2 = 0.5). In the nonthermalized 
chain (T = OK), the solitons attract one another, and the col- 
lision leads to their recombination. In the thermalized cyclic 
chain (T = 300K), their collision always leads to reflection 
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FIG. 16: Reflection of topological soli tons with different charges and 
polarities (qi = (1, 1);, q2 = ( — 1, l) r , si = — «2 = 0.5) in the 
thermalized (T = 300K) cyclic homogeneous AT chain (K = 0.234 
N/m, t r = oo). 



(Fig. I16> . This behavior can be explained by compression of 



phonons gas between solitons when they are drawing together. 
The compression leads to the repulsion of solitons, which in- 
creases as far as they are drawing together. In the chain with 
free ends, the compression of the phonon gas leads to long 
range repulsion of the solitons from the ends of the chain. 

Thus, topological solitons of the DNA chains are stable 
with respect to thermal oscillations. Interaction with thermal 
phonons does not lead to destruction or to the braking of the 
soliton, it leads only to changing the interaction between the 
solitons. In the thermalized chain, long range repulsion be- 
tween solitons is appeared. 



IX. CONCLUSION 

Investigation carried out in this paper shows that three types 
of topological solitons which imitate localized states with 
open base pairs, can exist in the considered asymmetrical 
model of the DNA double chain. It was shown that the soli- 
tons can move along the macromolecule with constant veloc- 
ity which is smaller than the sound velocity. In the inhomo- 
geneous chain, the character of the soliton movement depends 
on the sequence of base pairs in the molecule. In the chain 
with random inhomogeneous sequence, solitons can move at 
a distance no more than several hundreds of base pairs. The 
results of numerical investigations show that the solitons are 
stable with respect to thermal oscillations. Interaction of the 
solitons with thermal phonons of the macromolecule does not 
lead to destruction or to the braking of the solitons. And only 
the character of their interactions changes. The drawing of the 
solitons together leads to their repulsion, which is explained 
by compression of phonon gas between them. 

All these results point out that topological solitons of this 
type can be used to explain the long range effects in the DNA 
macromolecule. 
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